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Abstract 

A complete analysis of proton recoil asymmetry in neutron decay in the first order of radiative 
and recoil corrections is presented. The possible contributions from new physics are calculated in 
terms of low energy coupling constants, and the sensitivity of the measured asymmetry to models 
beyond the Standard model are discussed. 
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I. INTRODUCTION 



The free neutron decay, being one of the simplest semi-leptotic hadron decay processes, 
is very important in the search for possible manifestations of new physics. The main ad- 
vantage of neutron decay is the possibility to describe the process with minimal theoretical 
uncertainties and, as a consequence, the possibility to interpret unambiguously experimental 
results. The set of experiments for measurements of the neutron lifetime and neutron decay 
correlations can be used to determine the weak vector coupling constant, to test the uni- 
versal 



ity of the weak interaction, and to search for nonstandard couplings (see, for example, 
and references therein). The detailed analysis of required experimental 
accuracy and sensitivity to new physics of different observables for standard setups in neu- 
tron decay experiments have been done in the paper 9|. However, the recent measurement 
[lo| and the new proposal to measure 11 1 the integrated asymmetry of recoiled protons in 
relation to the direction of neutron spin (which is known as a C-angular correlation coef- 



ficient 



131]) raise the question about the sensitivity of this asymmetry to new physics. 



In order to be able to estimate the potential sensitivity of the C asymmetry to new physics 
and the best accuracy of the measurement of Standard Model parameters (e.g., the ratio of 
axial- vector and vector coupling constants of weak interaction), one needs to calculate recoil 
and radiative corrections for the C-asymmetry, as well as all possible contributions from the 
model beyond the standard one. Moreover, all these calculations must be done in the same 
framework to keep all possible uncertainties under control. 

In this paper, we use results of the effective field theory description of neutron beta- 
decay [l^ as a framework for the calculation of the C correlation coefficient in the Standard 
Model (with recoil and radiative corrections). Then we calculate possible corrections from 
new physics using the most general non-standard beta-decay interactions. This provides 
a consistent description of the proton recoil asymmetry in terms of low energy coupling 
constants related to models beyond the Standard one at a level well below that anticipated 
in the next generation of neutron decay experiments. 
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II. PROTON ASYMMETRY IN THE STANDARD MODEL 



We have chosen results, based on the effective field theory (EFT) approach, of the de- 
scription of the polarized neutron decay since this approach provides a general expression 
for neutron decay distribution function with the accuracy of 10~^ in terms of one free pa- 



rameter - low energy constant (LEG) (for more details, see paper To calculate the 

angular correlation coefficient C with the complete set of recoil and radiative corrections, 
we use a general expression for the differential neutron decay rate given by Eq.(8) in fl^ . 
It should be mentioned, that in the tree approximation (neglecting recoil corrections and 
radiative corrections), the EFT results reproduce exactly well known formula for neutron 



decay rate 



15| in terms of the angular correlations coefficients a. A, and B: 



^E,)Gl\Vud\\l + SX^' 



x(l + 6— + a ^ + A — — + B — (1) 
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Here, o is the neutron spin; is the electron mass, i?e, Ey^ and py are the energies 
and momenta of the electron and antineutrino, respectively; and Gp ys, the Fermi constant 
of the weak interaction (obtained from the /i-decay rate). The function $(-E'e) includes 
normalization constants, phase-space factors, and standard Coulomb corrections. For the 
Standard model the angular coefficients depend only on one parameter A = — C^/Cy > 0, 
the ratio of axial-vector to vector nucleon coupling constant (in general, Cy = C'y and 
Ca = C'j^ are low energy coupling constants for the low energy effective Hamiltonian given 

by Eq.my- 

1 - A2 ^ A2 - A „ A2 + A 

a= A = -2 — , B = 2 — . (2) 
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(The parameter b is equal to zero for vector - axial- vector weak interactions.) 



The C angular coefficient (do not mix with Cy and Ca) has been defined 12|] as the 
angular distribution of the recoil protons in the relation to the direction of the neutron spin, 
provided all other variables, including proton recoil momentum, are averaged out. In the 



tree approximation, it has been calculated in papers |12|. Il3[| . and numerical corrections to 
this approximation have been calculated in the paperjlj. Using this definition, one can 
calculate C-coefficient from a general expression for the differential neutron decay rate (Eqs. 
(8)-(19) in paper |l^) with all (in the first order) recoil and radiative corrections. To do 
this, we use the momentum conservation condition Pu + Pe + Pp = ^ which is multiplied by 



the neutron spin results in 



IPi^I cosOy + |pe| COS + IppI cos 6'p = 0, (3) 

where Pp is the proton momentum, and 9^, 6e, Op are angles between neutron spin and 
directions of anti-neutrino, electron and proton momenta, correspondingly. From Eq. ([3]) one 
can see that protons are going to the upper hemisphere (cos 6'p > 0),if cos6'i^+|pe| cos6'e < 
0, and to the lower hemisphere (cos^^p < 0), if Ip^l cos9u + \pe \ cos 6*6 > 0. Therefore, the C- 
coefficient, being a normalized difference of the neutron decay rate integrated over neutrino 
and electron angles, must be integrated over the electron energy under these two conditions. 
The integration over azimuthal angles leads to the Att"^ factors. To calculate integrals over 
6^ and 6^, it is convenient to work in cos-variables: co^Oy and cos^e- Thus, these two 
integrals could be represented in terms of a two-dimensional integral in (cos6'j^, cos^^g) space, 
which must be taken separately over lower and upper parts of the square area in the cosine 
plane: ([-1,1], [-1,1]). The line, dividing the area in two parts, is given by the equation: 
\py \ cosOy + \pe \ cosOe = 0. It should be noted that for both these integrals there are two 
different regimes of integration: > |pe| and \py\ < \pe\. For the first case, the integrals 
should be taken first over cos 6*^ and then over cos^^e, and for the second one in the opposite 
order. Applying this procedure for the decay rate given by Eq.([T]), one obtains 

where 



Xi = 5((E,"^"^)^) - 6(E™'^^)^m^ + 3E,""^m^ + 2m°/E™'^^ + 12E™'^^m^ ln(£;™'^^/me)). 

Eq.(jll) exactly reproduces the results of calculations of papers [l^, Q (the coefficient 
in Eq.(jl]) has a different sign since we define the positive direction of recoil protons as the 
direction of the neutron spin polarization). To obtain a general expression with radiative 
and recoil corrections, one has to apply the same procedure for the general neutron decay 



rate given by Eqs. (8)-(19) in paper 1^. These calculations are rather cumbersome but can 



be done exactly, without any approximation. Then, one can represent all corrections to the 
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C-coefficient in Eq.(jl]) as a sum of three terms 

AC = AC„ + ACs + ACrec, (5) 

where ACq contains Coulomb and radiative corrections, which do not depend on the nucleon 
structure (they are also known as the "outer" corrections ), ACs is the part of radiative 
corrections that is dependent on the nucleon structure (or the "inner" corrections), and 
ACrec represents recoil corrections. For recoil corrections we have 
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ACr 



12m„X(l + 3A2) 

X {9A/iy[(Er")' - m,')(Er")^ - 4(Er")'m2 + Smt + Amt ln(Er"Me)] 
+ A[31(E™"")^) - 117(E™"")^m2 + 279(^r")'^e - 211m^ + 18m^/(Er")' 
- 12m^(9(Er")' + llrnl) ln(Er7me)] 
+ 3/iy[-(E™'^^)^m2 - 9{E^^^ fmt + 9ml + ml/{E^^'')^ 



pmax-^e^ + 17{E;r^^yml + 9{E'^''''fml - 25m^ + ml/ {E'^^^'f 



- Um'MErr + ml) HErVm,)] (6) 
6XiA 



E™"V(E™«^)2 - ml[X\52{E^^^y - 12A{E^''^Yml + 507m^) 



5X(1 + 3A2 
+ (12(Ef "^)^ - 4{E^''''fml + 277ml)]}, 

and for strong interaction dependent part of the radiative corrections 

~ "2X(1 + 3A2)2^"^' 

where Cy is low energy constant (LEG) of the EFT [l^. The expression for ACa is very 
long and complicated to be presented here. However, one observes that all coefficients in the 
expressions for AC (including ACa) depend only on the mass of electron and the maximal 
electron energy. Therefore, one can re-write these expressions in a simple form (and without 
a lost of accuracy) by replacing the mass of electron and the maximal electron energy with 
their values: me = 0.511099 MeV and E™"^ = 1.293332 MeV. Then all dependencies on 
these parameters collapse to numerical coefficients in the front of neutron decay variables 
and the complete set of corrections AC could be written as: 

AC = -— -[;^(23.19375A + 4.45619A2) + ^e(5(0.2748 - 1.0993A) (8) 
(1 + 3A2) 271 271 



+ — (2.25672A - 0.265737A2 - 0.0113986/xv - 0.583714A/iy) 

rrin 



^ -(3.1326 + 7. 775A2)], 



m„(l + 3A2) 

where neutron mass m„ is in MeV. The first term in the first hne of the Eq.Q is ACq,, the 
second term is ACs, and last two hnes are recoil corrections. Now, using m„ = 939.57 MeV, 
fj,v = 3.7, a = 1/137.036, and A = 1.2695, one obtains 

AC = 0.0065 - 0.00022e(5. (9) 

Thus, all radiative and recoil corrections are expressed in terms of only one unknown pa- 
rameter (the EFT low energy constant) - which is supposed to be obtained from another 
independent experiment, if possible, or should be calculated from basic principles (for ex- 
ample, in lattice QCD). In the framework of the EFT, it could be estimated as Cy ~ 20 (see 
for details [l^). Discussions of another way of the estimation of Cy and its accuracy is given 
in the last section. 



III. NEUTRON /3-DECAY BEYOND THE STANDARD MODEL 

Now, when we understand all contributions to C angular correlation from the Standard 
model, we can consider how possible contributions from new physics can change the value 
of the C asymmetry. To calculate the possible contributions to the C-coefficient from the 
models beyond the Standard model, one can use the most general form of the Hamiltonian 
for the description of neutron /3-decay in terms of low energy coupling constants Ci (do not 
confuse with C angular correlation coefficient) by Q, 

Hint = {i'pi'n){Csi^ei^u + C'si^el5i'u) 

+ {i'p'l5i^n){Cp1pe'j5'^u + C'plpelpu) (10) 

+ Hermitian conjugate, 

where the index i = V , A, S , T and P corresponds to vector, axial-vector, scalar, tensor and 
pseudoscalar nucleon interactions. In this presentation, the constants Ci can be considered 



as effective constants of nucleon interactions with defined Lorentz structure, assuming tfiat 
all high energy degrees of freedom (for the Standard model and any given extension of 
the Standard model) are integrated out. Since we are interested in C angular correlation 
coefficient, which is the time reversal conserving one, all constants Cj can be chosen to be 
real. 

To see explicitly the influence of a non-standard interactions on the C angular coefficient, 
we will follow the procedure described in paper joj. First, we re-write the coupling constants 
Ci as a sum of a contribution from the standard model ^'^ and a possible contribution 
from new physics dCi. 



Cv 


— 


+ 5Cv 


C'y 


riSM 

— 


+ 5C'y 


Ca 


— 


+ 5Ca 


C'a 


riSM 

— '^A 




Cs 


= SCs 




C's 


= SC's 




Ct 


= SCt 






= 5C'j'. 





(11) 

I — I 

The pseudoscalar coupling constants are neglected here, since we treat [15|] nucleons non- 
relativistically. Then, we apply the described above procedure to the calculation of the 
C angular correlation coefficient from the Hamiltonian (ITOl) using Eq.(5) of paper Q]. (It 
should be noted that in the case of all 5Ci being equal to zero, the results is Eq. (jlj).) The 
obtained corrections to the C correlation coefficient due to contributions from non-standard 
modes can be written as: 



SC^ewPhys — TTTTTZ — . „ , ox + 



2X(1 + 3A2) 2X(1 + 3A2) 



Xi2\ 



X(1 + 3A2) 

where 



+ 



:i + 3A2) X(1 + 3A2) 



(12) 



X2 = 10me(Er") V (^I""")' - m2(2(Er")' + 13m2) 

- 30me(4(Er")'m,' + mt) In ((E^" + sgrtiE^"? - ml)/me) 
X3 = 5me(3(E™"^)^ + 12{E^''''fml - 15m^ - 12{2{E^'''=fml + m^) ln{E^^yme)). 
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The coefficients Li depend only on new physics contributions: 



{6Cv + 5C'y) + {5Cv^ + 5C'y^ + SCs" + 5C^')/2 



L 



+ ?>[\{5Ca + 5C'j^) + {5Ca^ + 5C'a^ + 5Ct^ + 5C^')/2], 
1 = -2{5Ca + 5C'a) + ?>5Ct5C't - 2[bCybC'j^ + + 2\[bCy + 



(13) 



(14) 



= yr^^[-2A(5Cr + 5C^) - K^Cs + 5C;;) + (5Ct + C^) 



(15) 



In the above expressions, we have neglected radiative corrections and recoil effects for the 
new physics contributions, but kept Coulomb corrections since they can be important for a 
low energy part of the electron spectrum. 

From Eq.( fT2l) . one can see that, as in the case of radiative and recoil corrections, all coef- 
ficients in the expression are functions only of electron mass and maximum electron energy. 
Therefore, we simplify the general expressions for the contributions from new physics, by 
substituting numerical values for all known parameters (electron mass, electron maximal 
energy, as well as for a = 1/137.036 and A = 1.2695) and keep only first order contributions 
from non-standard interactions. Then, Eq. f|T2l) transforms into 



Instead of the presentation of these corrections in terms of low energy coupling constants 
related to the Lorentz structure of weak interactions, we can re-write them in terms of quark 



SCNewPhys = 0.05657(5Cy + 5C'y) + 0.04456(5^^ + 6C'^) 
- 0M234{6Cs + 6Cs) + 0.02132{6Ct + 5C't). 



(16) 



and lepton current constants aji and Aji^ defined in paper 




[3] . Using the transformation rules 



5Cy + 5C[ 



V — 



2{aLL + CtLi?), 



5Ca + 6C'a = 2\{aLL-aLR), 
6Cs + 6C's = 2gs{ALL + ALR), 
6Ct + SC^ = Aqtcill, 



(17) 
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TABLE I: Possible manifestations of new physics 



Model 


L-R 


Exotic Fermion 


Leptoquark 


Contact interactions 


SUSY 


Hiffgs 


a-LL 






0.2 - 0.03 








aLR 




0.01 


0.01 








All + Alr 








0.01 


7.5 • 10-4 


3 • 10-6 


-All + Alr 






3 • 10-6 









and assuming j?! gs 
physics as: 



1 and = 1, we obtain the expression for corrections from new 



5C NewPhys = 0.11314(aLL + a^/j) + 0.11314(aLL - cllr) 
- 0.12A68iALL + ALR) + 0.08528aLL- 



(18) 



The parameters dji, aji, and Aji describe contributions to the low energy Hamiltonian 
from current-current interactions in terms of j-type of leptonic current and i-type of quark 
current. For example, cilr is the contribution to the Hamiltonian from left-handed leptonic 
current and right-handed quark current normalized by the size of the Standard Model (left- 
left current) interactions, gs and gT are formfactors at zero-momentum transfer in the 
nucleon matrix element of scalar and tensor currents. For more details, see paper Q]. 

The expected values of these parameters vary over a wide range from 0.07 to 10-6 (see 
Tableland paper j?] for the comprehensive analysis and for discussions of significance of 
each of these parameters for models beyond the Standard one). 



IV. CONCLUSIONS 

Taking into account the results of Eqs.(jl]), ([5]), f|T2|) . and (|T8l) . one can write the complete 
expression for the C angular coefficient (Ctotai) as a sum of the tree-level approximation C, 
radiative and recoil corrections in the Standard Model AC, and possible contributions from 

new physics 6C NewPhys- 

Ctotal = C + SC + 5C NewPhys- (19) 

It should be noted, that this equation is the exact expression of the C angular correlation 
coefficient in the first order of recoil corrections, radiative corrections, and low energy con- 
tributions from new physics. Therefore, it could be considered as the complete expression 
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up to the level of accuracy of 10~^, provided the EFT low energy constant (LEG) is given. 
Otherwise, it could be considered as a parametrization in terms of one free parameter - 
LEG with the same accuracy of 10"^. Would the parameter Cy be determined from another 
independent experiment (for example, from the precise measurement of neutrino-deuteron 
cross-sections) or calculated using lattice QGD approach, Eq. (HM could be used to test 
the Standard model up to the level of accuracy of about 10~^, by comparing a theoretical 
prediction with experimental results. Unfortunately, neutrino experiments and QGD calcu- 
lations with the required accuracy are rather difficult problems and we cannot rely on them 
at the present time. 

To understand the desirable level of accuracy in a search for new physics, one can use 
first a conservative approach: the estimate for the LEG as Cy ~ 20 given in paper fl^ . 
Then, the level of theoretical uncertainties due to strong interactions, according to Eq.Q, 



is about 0.0044, w: 



fich is comparable to the claimed experimental accuracy 0.0026 of the 



recent experiment 10|. However, as it was mentioned in iJ] that by comparing the results 



of the EFT approach and the calculations of radiative corrections for total neutron decay 

n n Q 

rate U7, Qj, Qj, one can find the correspondence between these two calculations, which 

3 



results 



in the following equation 



.^-4inr!^U3inr!^Uinr^ 

4 \mz / \mj^J KniA 



2C 



Born 



+ Aa. 



(20) 



Here mwi^z are the masses of the W, Z bosons and tua is the axial mass scale, which 



are rather well known. The source o 
CBorn and Ag (see, for details 



17 



18 



theoretical uncertainties is related to two last terms 



19|, 



20| ) . Ghanging from the EFT " ideology" with one 



unknown LEG to direct calculations using strong interaction models, we lost the attractive 
feature of the model independent EFT approach and have to deal with dependencies on 
strong interaction models applied for description internal structure of nucleons. On the 
other hand, in the given framework 17, [is], [3]) which is actually a very well recognized 



standard approach to general analysis of weak interactions, we can reduce uncertainties in 
the estimation of LEG to the uncertainties of calculations of CBom and Ag terms. Then, 
using results of recent calculations of these terms |20|] CBom — 0.829 and Ag ~ —0.34 
with the claimed level of uncertainty of 10%, one can reduce the level of uncertainty of the 



obtained theoretical description of the C angular coefficient to the level of about 10 



1. e. 



to the level of validity of the description of neutron decay in paper 
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Accepting these estimates, one can see from Eqs. (fT6|) and (fTSi) that precise measurements 
of the C angular correlation can provide limits for non-standard interactions in terms of 5Ci 
coupling constants up to the level of about (2 — 5) • 10~^, or, in terms of parameters related 
to non-standard currents, up to the level of about 10^^. However, in order to be able to 
constrain new physics parameters at this level, the currently achieved experimental accuracy 
jlo| must be improved by two orders of magnitude. 
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